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Abstract 

We discuss a class of time-dependent Hamilton-Jacobi equations, where an unknown function of 
time is intended to keep the maximum of the solution to the constant value 0. Our main result is 
that the full problem has a unique viscosity solution, which is in fact classical. The motivation is a 
selection-mutation model which, in the limit of small diffusion, exhibits concentration on the zero 
level set of the solution of the Hamilton-Jacobi equation. 

Uniqueness is obtained by noticing that, as a consequence of the dynamic programming principle, 
the solution of the Hamilton-Jacobi equation is classical. It is then possible to write an ODE for 
the maximum of the solution, and treat the full problem as a nonstandard Cauchy problem. 
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1 Introduction 

1.1 Model and question 

The purpose of this paper is to discuss existence and uniqueness for the following problem, with 
unknowns x))\ 

dtu = |Vup + R{x, I) {t > 0,x G maxu(t, x) = 0 

X 

< 7(0) = Jo > 0, (1) 

u{0,x) = uoix), 

where Iq > 0 and uq is a concave, quadratic function. For a given continuous function I{t), u{t,x) 
solves a Hamilton-Jacobi equation. The unknown I may be thought of as a sort of regulator, or a sort 
of Lagrange multiplier, to maintain the maximum of u equal to 0. The constraint on the maximum 
of u{t ,.) makes the problem nonstandard. 
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Existence of a solution [u, I) to ([T|) is not new, the first result is due to Perthame and Barles [2T] (see 
also Barles-Perthame-Mirrahimi [3] for a result with weaker assumptions). An important improvement 
is given by Lorz, Perthame and the first author in [16]; they indeed notice that a concavity assumptions 
on R - that we also make here - entail regularity. This allows them to derive the dynamics of the 
maximum point of a solution u{t,x). See also [T7]. Both types of results rely on a special viscous 
approximation of m - see equation (|4|) below. Uniqueness, however, has remained an open problem, 
apart from a very particular case m- 

The main goal of the paper is to prove the missing uniqueness property; a result that we had already 
announced in |19] . We also provide a constructive existence proof which was not available in the pre¬ 
vious existence results Two important consequences, that we will present in a forthcoming 

paper [20] (see also |19]). will be the convergence of the underlying selection-mutation model in a 
stronger sense than what is known, and asymptotic expansion of the viscous solution. The asymptotic 
expansion, which allows to approximate the phenotypical distribution of the population when the mu¬ 
tation steps are small but nonzero, is particularly interesting in view of biological applications. One 
of the main ingredients will be regularity under suitable concavity assumptions on R and uq, which is 
far from being available in general. Instead of relying on viscous approximations we will prove these 
results directly for the equation 

ut = \Vuf + R{t,x). (2) 

This will allow a much easier treatment than in the usual viscosity sense. 

The uniqueness result will also be helpful to develop the so-called Hamilton-Jacobi approach (see for 
instance Diekmann et al. m, mm and Subsection II. 2h to study more complex models describing 
selection and mutations. For instance, our result would allow to generalize a result due to Perthame 
and the first author in |18] on a selection model with spatial structure, where the proof relies on the 
uniqueness of the solution to a corresponding Hamilton-Jacobi equation with constraint. 

1.2 Motivation 

Model ([T]) arises in the limit e —>• 0 of the solutions to the problem 

dtUe - eAris = —R{x,Ie{t)) (t > 0, x E R*^), /^(t) = [ ii{x)ne{t,x)dx, (3) 

where n^{t,x) is the density of a population characterized by a d-dimensional biological trait x. The 
population competes for a single resource, this is represented by leit), where tp is a given positive 
smooth function. The term R{x, I) is the reproduction rate; it is, as can be expected, very negative 
for large x and decreases as the competition increases. The Laplace term corresponds to the muta¬ 
tions. The small parameter e is introduced to consider the long time dynamics of the population when 
the mutation steps are small. Such models can be derived from individual based stochastic processes 
in the limit of large populations (see Champagnat-Ferriere-Meleard mm- There is a large litera¬ 
ture on the models of population dynamics under selection and mutations. We refer the interested 
reader, for instance to Geritz et al. m and Diekmann [8| for an approach based on the study of 
the stability of differential systems (the so-called adaptive dynamics approach), to Champagnat |1| 
for the study of stochastic individual based models, to Raoul [22] and Mirrahimi m for the study of 
integro-differnetial models. 

The Hopf-Cole transformation = exp (us/s) yields the equation 

dtUe = sAUe + lVUel"^ + R{x,Ie) (4) 
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which, in the limit e ^ 0, yields the equation for u. Furthermore, Ig being uniformly positive 
and bounded in e, the Hopf-Cole transformation leads to the constraint on u. One expects that 
concentrates at the points where u is close to 0 and the function appears, in the limit, as a sort of 
Lagrange multiplier. One has indeed 

ne{x,t) —^ n{x,t) = p{t)6{x — x{t)), weakly in the sense of measures, 

with 

u(t,x(t)) = maxuft,x) = 0, pit) = S}.. ■ 

X 'ip{x[t)) 

This method to study ([3]) has been introduced in [9] and then developed in different contexts. See for 
instance Perthame-Barles [21] (convergence to Hamilton-Jacobi dynamics for Q, Barles-Mirrahimi- 
Perthame 0 (the same type of result, but with nonlinear, nonlocal diffusion), Champagnat-Jabin 
[7] (nonlinear integro-differential model with several resources), |16] (convergence improvement by 
introduction of the concavity assumptions). This approach has a lot to do with the ’approximation 
of geometric optics’ for reaction-diffusion equations of the Fisher-KPP type. See Freidlin [111 |T2] for 
the probabilistic approach, and Evans-Souganidis [TO], Barles-Evans-Souganidis [2| for the viscosity 
solutions approach. 


1.3 Assumptions 

The assumptions we are stating below are in the same spirit (but slightly weaker) as in [16], where the 
authors noticed that this set of assumptions allowed them to work with smooth solutions, thus going 
quite far in the study of ([3|). We believe that the results that we will prove would certainly be false if 
some of those assumptions were removed. 

• Assumptions on R{x,I). We choose R to be smooth, and we suppose that there is Im > 0 
such that (fixing the origin in x appropriately) 


maxi2(x, Im) = 0 = 12(0, Im), 


(5) 


“ < -^(3^) 7) < Aq — iFi|x|^, for 0<I<Im, 

— 2K] < D^R{x,I) < —2Ki < 0 as symmetric matrices. 


( 6 ) 

(7) 

( 8 ) 


d^R d^R 

'(^^7)|-F| (x,/)|<A 3 , for 0 </< Im, and f, j = 1,2,-•• ,f2, (9) 


didxj 


didxidxj 


\D^Ri;I)h. 


< A4, 


for 0 < / < Im- 


( 10 ) 


Assumptions on uq{.) and Iq. We assume the existence of positive constants Lq, Lq, L^, Li 
such that 


Note that this implies 


~ Lo ~ Lli\x\^ < uo(x) < Lq — Ti|x|^, 
— 2T]^ < D^uq < — 2Ti. 

\Duq{x)\ < L2(1 -k |x|). 


( 11 ) 

( 12 ) 

(13) 
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for a large constant L 2 > 0. We also need that, for a positive constant L3, 


\D^ 




< L3. 


(14) 


Finally we assume that 


max uo{x) = uo{xo) = 0, R{xo,Io) = 0. 

X 


(15) 


Note that the monotony assumption ([5]) means that the growth rate decreases as the competition in¬ 
creases, which is natural from the modeling point of view. The concavity assumption is a technical one. 

In Section [2] we will study an unconstrained Hamilton-Jacobi equation where we replace R{x,I) by 
R{t,x). To prove our results on this unconstrained problem we assume same type of regularity and 
concavity assumptions on R that we state below: 

• Assumptions on R{t,x). We choose R to be smooth, and we suppose that 

— K.i\xf ^ < Kq — Ki\x\‘^, for t E M^, 

— 2K I < D‘^R{t,x) < —2Ki < 0 as symmetric matrices, 

\\D^R{t, •)llL°°(Rrf) ^ -^4) ^ 

1.4 Results and plan of the paper 

Our first result concerns the unconstrained Hamilton-Jacobi equation 

ut = I Vup -|- R{t, x) (t > 0, X E M"^), 
u(0, x) = uo{x). 

The assumptions on uq and R are those stated in the preceding subsection. 

Theorem 1.1 Equation (|19p has a unique viscosity solution u that is bounded from above. Moreover, 
it is a classical solution: u E lFioc°°(IR'^)) nlFj]j’“(M'''; Lj^^(M'^)), — max(2L;^, \/K]) < D‘^u < 

— min(2Li, -s/Ki) and ||T>^'u||looqo,t]xr 4) < Li{T) where L/i{T) is a positive constant depending on 
Lli^ K.I 1 -^4) -^3 T. 

Let us point out that the assumption of boundedness from above is, most certainly, irrelevant. However 
the constraint in ([1]) implies that all the solutions that we consider are bounded from above. This 
extra condition is thus legitimate, and will keep the length of the preliminary work to a minimum. 

Theorem 1.2 The Hamilton-,Jacobi equation with constraint 0 has a unique solution {u,I). More¬ 
over we have 

{u,I) E L^,{R+-W^^{R<^)) nW^^{R+-L^,{R'^)) x ^^’-(M). 


(16) 

(17) 

(18) 


(19) 


The paper is organised as follows. In Section 2 we prove the Cauchy Problem for (jl9p . In Section 
3 we reduce ({IJ to a (nonstandard) differential system. Theorem 11.21 is proved in Section 4. Section 5 
is devoted to the study of a particular example. 
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2 The Cauchy problem 

In this section, we prove Theoreni ll.il 

We will first prove that the only solution to ([2]) that is bounded from above is the solution u{t, x) of 
the dynamic programming principle 

u{t,x)= sup \f{j) : 7 E , (20) 

(7(s),s)eK'ix[0,t] 

with 

F{j) := no(7(0)) + ^ + R{s,j{s))^ ds. 

Such a solution is a viscosity solution to ([II). We will prove, in addition, that it is classical and satishes 
the properties claimed by Theorem ll.il 


Uniqueness. This step essentially consists in showing that a viscosity solution of ([2]) does not grow 
too wildly, which will reduce the problem to the application of classical arguments. We have already 
assumed boundedness from above, so let us show that u goes to —oo at most in a quadratic fashion. 
Due to the assumptions (fT6]l on R, we have 

dtv>0, v(t. x) = u(t. x) + tK^\x\^ 


in the viscosity sense, which implies that v is time-increasing - thus the needed estimate. Let us - 
although this is elementary - explain why: choose T > 0 and assume the existence ofO<s<t<T 
such that the inequality v{t, x) > v{s, x) does not hold. In other words there is xq such that v(s, xq) > 
v{t,xo). For e > 0 consider the quantity 


We{t,x) 


v{t,x) - v{s,xo) -I- 


T-t 




\x - Xol^ 


For e > 0 small enough, there is a local minimum point for Ws, called Xe), such that te is bounded 
away from s and T, and —>• xq as e —>■ 0. At that minimum point the viscosity inequality implies 



> 0 , 


a contradiction. 


So there is a large C > 0 such that 


—(7(1 -|- t)(l -|- |xp) < ti(t, x) < C. 


And so, by an easy adaptation of Chap. 2 of Barles [T], where a uniqueness result for solutions that 
grow at most exponentially fast is provided (see also mm), there is at most one viscosity solution 
to (l 2 |) that is bounded from above. 


Existence. We may thus turn to (1201) . Let us suppose that (7n)i<n) with 7 „ E (7^([0,t];M'^) and 
7n(t) = X, is such that F{'yn) —t u{t, x) as n —oo. Since R and uq are bounded from above, we obtain 
that, for some constant C 

\jn\‘^{s)ds < C. 
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Consequently, from 7 n(t) = x we deduce, modifying the constant C if necessary, that 


ll7n||l^l,2[o,t] < C. 


It follows that, there exists 7 E Vh^’^([0,t];M'^), such that as n — 00 , 7 ^, —)• 7 strongly in C([0,t];M'^) 
and weakly in Vh^’^([0, t]; M'^). We deduce that, as n ^ 00 , 


Wo(7n(0)) Mo(7), / R{s,-in{s))ds 


f 


R{s,-f{s))ds, 


f 


I7P 


{s)ds < liminf 

n^oo 


|7n|^(s)ds. 


We conclude that ^ ^ 

u{t,x)= ^ 0 ( 7 ( 0 )) + ^ + i 2 (s, 7 (s))^ ds. ( 21 ) 

We claim that such a trajectory is unique, which entails uniqueness for the Cauchy problem (fT9]l . We 
note indeed that such trajectory 7 satisfies the following Euler-Lagrange equation 


' 7 ( 5 ) = - 2 Vi?(s, 7 (s)), 

< 7(0) = -2Vuo(7(0)), ( 22 ) 

j{t) = X. 


However, from the concavity assumptions on R and uq we obtain that the above elliptic problem is 
coercive and hence the solution 7 is unique. 


Regularity. Let us denote by ^x{t) the unique solution of (l^2|) . The function (t,x) e-)- 
belongs to L5^^(M+, ^^^^(M'^)) because Vuq and VR{t,.) are in Now, we have 

u{t, x) = no(7a;(0)) + E(7a;), 

yielding (this is a classical computation): 

diu{t,x) = Vrio(7x(0)).9i7x(0) J (^~ R{s,jx{s))dijx{s)^ ds. 

Integrating by parts and using the Euler-Lagrange equation ([2^ . we get: 

Vuit,x) = -^. (23) 

Thus u E (M+, 

Strict concavity. We will prove that u{t, x) is uniformly strictly concave, namely that D^u < —2X1 

in the sense of symmetric matrices, for A = min (Li, ^'^^ ). To this end, we show that, for all a E [0,1] 
and (x, y) E X 

au{t, x) -I- (1 - a)u{t, y) + A(t(1 - a)\x - y|^ < u{t, ux + (1 - a)y). (24) 

Let 7 a; and jy be optimal trajectories, solving (1221) . with 'jxit) = x and 7 j/(t) = y. Note from the 
choice of jx and 7 ^ that we have 

u{t,x) = uoi'jxiO)) + j + R{s,^x{s)^ ds, 
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and 


i{t, y) = uo( 7 j/( 0 )) + ^ + R(^s, 7 y(s))^ ds, 


u{t, ax + {l- a)y) > uo{a-fx{0) + (1 - o-) 7 ?;( 0 )) 

f k7x + (1 - o')%(s)P 


/ 


+ R{s, ajx{s) + (1 - o-) 7 y(s))^ ds. 


Furthermore, from the concavity assumptions on R and uq we have 

auo{t,x) + (1 - a)uo{t,y) + Lia{l - o-)| 73 ,( 0 ) - 7 j/( 0 )p < uo{t,ax + (1 - a)y), 

and ^ ^ ^ 

a / i?(s, 7 a;(s))ds + (1 - fj) / R{s,-fy{s))ds + Kia{l - a) / 173 ,( 5 ) - 7 y(s)pds 
Jo ^ Jo Jo 

< / R{s,ajx{s) + (1 - o-)jy{s))ds. 

Jo 

Moreover, from the strict concavity of we obtain that 


■i 


a / - 


‘ l7.(s)|^„ , ,3 _, /■* l7.(s)P,„ , /■M7.(^)-%(^)P 


4 4 


-ds + a(l — 


I 


< - 


\(T% + {I - a)%{s)Y 


{l-a) f 

Jo 


ds 


ds. 


We deduce that 


u{t, fix + (1 - a)y)) > au{t, x) + (1 - cr)u{t, y) 

+ - c') f / (7l7x(s) - %(s)P + di:i|73,(s) - 'yy{s)\^)ds + Fi|73,(0) - jy{0)Y 


Next we have 

V% d_ 
2 Jo ds 

Writing 


l7x(s)-7y(s)Pds < |73 ;(s)- 7y(s)pds + ^ i^^^-^^(s)ds. 


\x - yP = hxit) - 7y(i)P = l7x(0) - ly{Q)\^ + - ly{t)fds 


(25) 


we find 

k-yP < l 7 x( 0 )- 7 y( 0 )P+ 2 y^^ |73,(s) -7j^(s)pds +—I73, - 7yP(s)ds. 
Combing the above line with ()25p . we obtain (1241) for A = min (Li, ■ 

Bounds on u and Vn. The first thing to notice is a bound for ( 73 ,, 73 ,)- Indeed, the coercivity for 
(1221) . as well as the fact that x 1 —>■ Vi?(t, x) grows linearly with a constant only depending on t, implies 
the existence of a locally bounded constant K(t) such that 


|(7x('S),7a:(s))| < Kit){l + la;|), for all s G [0,t]. 


(26) 
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This implies, modifying the constant K{t) if necessary, thanks to (I2ip : 

\u{t,x)\ < K{t){l + |x|^). 


Moreover, from (|2c{l) . we have 


\Vu{t,x)\ < K{t){l + |x|). 


(27) 


Semi-convexity. Because u is three times differentiable in x, then dtu is locally in x and 

equation (fT9]l may be differentiated twice with respect to x. So, let e be any unit vector, we have 


dtideeU) = 2|V5eUp + 2Vu.V(5eeU) +deeR{t,x) 
> 2\deeU\^ + 2Vu.V(9eeU) + deeR{t, x). 


Because of (l27)l . the curve 1'Jxit) becomes a characteristic curve for the equation 

dtv = 2v^ + 2Vu(t, x).Vv(t, x) + deeR{t, x). 


Moreover, along this characteristics, we find 

= 2v{t,-fxit)f + deeRit,'yxit)). 

We deduce, thanks to (fT^ and (flTll . that 

deeU > — max(2 L^, 

Bounds on D^u. As for the third derivative, we set v{t,x) = D^u{t,x), the equation for v - that we 
may obtain using differential quotients - is 

dtv — 2Vu.Vu = S{t, X, v) := Qv.D^u + D^R, 

where Vu denotes the column of tensors {div, ...,ddv) and v.D^u denotes the column of matrices 
{diD'^u.D‘^u, ...,ddD^u.D^u). This is a linear equation with (thanks to the bound on D‘^u) bounded 
coefficients. Thus, local boundedness of ||u(t, ■)||ioo(]gd) holds, and this concludes the proof of Theorem 

o 


3 Uniqueness: reduction to a differential system and properties of 
the function I 

The idea is to change the constrained problem ([T|) by the following slightly nonstandard differential 
system: 

R{x{t),I{t)) = 0, for t E R+, 

< x{t) = [—D‘^u(t,x{t))) ^ Vi?(x(t),/(t)), for t € M+, 

^dtu = \Vu\‘^ + R{x,I), in R+ X 

with initial conditions 

7(0) = Jo, u(0, •) = iio(-)> x(0)=xo, 

such that m.aXxUo{x) = uo{xo) = 0 and R{xo,Io) = 0. 

Note that (1281) is really a differential system because the assumptions on R imply that I{t) can 
implicitely be expressed in terms of x{t). And it is slightly nonstandard because x solves an ODE 
whose nonlinearity depends on u. The precise statement is the following 


(28) 


(29) 



Theorem 3.1 Solving the constrained problem ([T]) is equivalent to solving the initial value ODE-PDE 
problem (f28]) - (f^ . 

Proof. Let (n, /) be a solution of ([T]) with initial datum {uq,Iq), the function I being continuous, 
and u a solution of the Hamilton-Jacobi equation in the sense of (120^ . Theorem o is applicable, and 
yields a solution u{t, x) which has at least three locally bounded spatial derivatives, locally uniformaly 
in time. Moreover, the D^u is bounded uniformly in time and in x, and finally the function u{t, .) is 
strictly concave. This allows a lot. 

• There is, at each time, a unique x{t) maximising u{t, .) over M'^. Thus the trivial identity 

Vu{t,x{t)) = 0 (30) 

can (use differential quotients) be differentiated with respect to t, to yield that (i) x{t) is locally 
W^’°° and (ii) the (a priori less trivial) identity 

dt{Vu){t,x{t)) + D‘^u{t,x(t)).^{t) = 0. (31) 

• The function u{t, x) has enough regularity so that we may take the gradient of ([I|) with respect 
to X, and evaluate the result at x = x{t). Because of (IHH]) we have D‘^u{t,x{t)).Vu{t,x{t)) = 0 
and, because of (l3T]l . we have 

- D^u{t,x{t)).^{t) = VR{x{t),I{t)). (32) 


The last item to take into account is the constraint 

u{t,x{t)) = 0, 

which we may (still with the use of differential quotients) differentiate with respect to time, in 
order to yield 

dtu{t,x{t)) + Vu{t,x{t)). — {t) = 0, 


dt 


thus entailing 
This yields, from ([1]), 


dtu{t, x{t)) = 0. 


R{x{t), I{t)) = 0. 

Gathering ([331) . ([I]) and (l32|) shows that the constrained problem implies 


(33) 


We also prove that regularity plus (l28|) easily implies ([T]). Let {u,I) solve ([281) . The first line of ([T]) 
derives immediately. To prove the second line, note that, thanks to Theorem 11.11 and the third line 
of (j28ll we deduce that u is strictly concave. It has hence, for all t E M+, a unique strict maximum 
point, in the variable x, that we denote y{t). Following similar arguments as above, we obtain that 

iy{i) = {-D'^u{t,y{t))~^ VR{y{t),I{t)), for t E M+, 

\y{0) = xq. 
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Comparing this with the second line of (1281) and (I29p we obtain that x{t) = y{t), for all t G M'*'. Finally, 
evaluating the third line of (|28l) at x{t) and using the first line of ([28]) . we obtain that dtu{t,x{t)) = 0. 
This equality together with Vu{t,x{t)) = 0 and (l2^ implies that 

maxu{t, x) = u{t,x{t)) = 0. 

X 

and the proof of Theorem 18.11 is complete. □ 

Remark. What we have done here is nothing else than the derivation of the equation for x carried out 
in m- In particular, the equilibrium property R{x, /) = 0 would hold in a more general setting than 
here. The new point here is that we establish, in a mathematically rigorous fashion, the equivalence 
between the initial problem dm and the coupled system (1281) . and this equivalence holds because of all 
the differentiations with respect to x and t that we are allowed to make. 


4 The proof of Theorem 11.2 

We fix T > 0. To prove that ([T|) has a unique solution {u, I) in [0, T] x it is enough to prove that 
there exists a unique solution to (I28l) ~ (f2^ . 

We prove this using the Banach hxed point Theorem in a small interval and then iterate. 


4.1 The mapping $ and its domain 

First, we define 

n = {x I R{x, 0) > 0} , rio = {x I R{x, Jo) > 0} , and 
A = |x(-) G C(^[0,(5];R(xo,r5)^ x(0) = xoj , 

where B{z,r) is the ball of radius r centered at z, and <5 is a positive constant such that 
(5 < min(/x, c(r)), with /x = — ' ^ ^ d{I}o , 11), 

Cm 

where d{A,B) is the distance between the sets A and B. The constant rg is given by 

CmS 


rs = 


min(2Li, s/l^) 


(34) 


the constant Cm is chosen such that 

|VR(x,/)| < Cm, inOx[0,/M], (35) 

and c(T), a constant depending only on T, will be chosen later. Note that, by the choice of 5 and 
since xq G flo, we obtain that 

B {xo,rs) C Q. 

Our theorem will be proved by the introduction of a mapping 

: .4 —)• .4, <h(x) = y. 


We will prove the following theorem. 
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Theorem 4.1 The mapping ^ is a strict contraction from A into itself. 


To define we first need to introduce some other mappings. Let x(-) G A. We define I : A ^ 
C ([0, (5]; [0,/m]) such that 

R{x{t),I[x]{t)) = 0. 

From ([5|), ([8|) and since x{t) G for all t G [0,(5], it follows that I can be defined in a unique way. 
Next, we define the following domain 

B= [ve L-([0,(5]; n 5]; | 

- max(2 Li, y^) < D‘^v < - min(2Li, \/^), II^^^^IIl°°([o,5]xR'^) < , 

and the following mapping 

fy :C([0,(5];[0 ,Im])^^ 

\V{I)=V, 

where v solves 

[dtv = \S/v\‘^ + R{xJ), in [0, (5] X M'^, 

|^u(0, x) = uo(x), in 

It follows from Theorem II.II that the above mapping is well-defined. 

Finally, we introduce a last mapping: 


IF: C([0,(5];[0,/m]) xAxB^A, 

\f{I,x,v) = y, 

where y G A solves 

iy{t) = {-D‘^v{t,x{t))~^VR{x{t),I{t)), in [0,(5], 

\y{0) = xq. 

To prove that F is well-defined, we must verify that y{t) remains in B {xo,rs). We note that, since 
V G B, we have 

0 < - 1 =^ < {—D'^v{t,x{t))) ^ ^- 

max(2 , y^) min(2Li, V^) 

We deduce, thanks to (l35|) . that 

y{t) G B {xo,rs ). 


We are now ready to define mapping <I>: 


$ : A —y A, 

$(x) = F(X(x),x,y(X(x))). 


It follows from the above arguments that $ is well-defined. 
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4.2 Proof of Theorem 11.21 


In this section we will explain why uniqueness to (|28p holds. One main technical lemma will be stated, 
its proof will be postponed to a special section. 

Let us prove that is a contraction for c(T) (and hence 6) small enough. To this end, we first prove 
that I is Lipschitz: 

|X(xi)-X(x 2 )| < C||xi - X 2 ||loc([o, 5 ]), for all xi, a ;2 G C(^[0, d]; S(xo, . (37) 

We have indeed 

R{xi,I{xi)) = R{x 2 ,T{x 2 )) = 0. 

It follows that 

R{xi,l{xi)) - R{x 2 ,I{xi)) = R{x 2 ,I{x 2 )) - R{x 2 ,I{xi)), 

and thus 

dR 

VR{cXl + (1 - c)X 2 ,/(xi)) • (X 2 - Xl) = -^{X 2 , J)(/(x 2 ) “ 

with c G (0,1) and J G (/(xi),I(x 2 )). Finally, using ([8|), ([35]) and the fact that xi(t), X 2 (t) G 0 for 
all t G [0, (5] we obtain (|3^ with C = 

Next, we have that V : C ([0, d]; [0, Im]) B is also Lipschitz. 

Lemma 4.2 Let R, R & C{[0,6];[0 ,Im])- Then 


||I/(/i) - P(/2 )||w2.<x>([o,5]xR‘*) R C\\Ii - I2 \\l^{[0,S])^- (38) 

This is a nontrivial lemma, whose proof will be given in the next section. 

From ([7|) and dH), and since xi, X 2 G .4 and W, V 2 G B, we deduce that F : C ([0,5]; [0, Im]) 'xAxB — 
A is Lipschitz with respect to all the variables with Lipschitz constant C6: 


||F(/i,Xi,Fl) - F(/2,X2,P2)||l<»([0,5]) 


< C6 


\xi - a;2||L°°([0,5]) + 11-^1 - -^2||l°°([0,5]) 


+ IIP1-P2 


2|liv2.°°([0,5]xll 


(39) 


Finally, we conclude from (f37|l . ([38l) and ([3^ that ^ 4. is a Lipschitz mapping with a Lipschitz 

constant C6: 

||4>(xi) - 4 >(x2)||l°°([0,5]) < C6\\xi - X 2 11 LOO ([0,5]). 

Choosing c(T) (and hence d) small enough, we deduce that is a contraction. 

We deduce from the Banach fixed point Theorem that 4* has a unique fixed point and consequently 
(l28l) - ([29]l has a unique solution for t G [0,5]. 


To prove that (I28l) - (l2^ has a unique solution in [0,r] we iterate the above procedure K = \j~\ 
times. Let 2 < i < K and (x,/, tt) be the unique solution of for t G [0, {i — 1)5]. Then, at 

the Tth step, we consider the same mapping 4> but as the initial data we choose 

xo = x((i - 1)5), So(-) = ^((*-1)'5>-), Io = I{{i-i-)S). (40) 

We claim that these initial conditions satisfy 


xo G Po) ^^0 £ — max(2L;^, -\/K)') < D'^uo < — min(2Li, \/1r) (41) 

llL>^no||L°°(Kd) < TRii - 1)5), ma.x^uo{x) = xq, R{xo,Io) = 0. 
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Lemma 4.3 Let {x,I,u) be the unique solution of (I28|) - (l29]) in [0,t]. Then, I{t) is increasing with 
respect to t in [0, r] and 

x(i) G Qq, M(t, •) G with — max(2 L^, ^K^i) < D‘^u{t, x) < — min(2Li, a/ Ki), 

||-D^u||£,ooqo,t]xR<*) ^ m.8LKxu{t,x) = xft), R{x{t), I (t)) = 0,, for all t ^ [0,t] andx^M!^. 

Proof. The regularity estimates on u are immediate from Section [2l The two last claims also follow 
immediately from (|28l) - (l29p and the arguments in Section [3l We only prove that I{t) is increasing 
with respect to t and x{t) G LIq. 

Differentiating the first line of (|28p with respect to t we obtain 

VR{x{t),I{t)) •T(t) + ^R{x{t),I{f))i{t) = 0. 

Moreover, multiplying the second line of (j28p by Vi? we obtain 

VR{x{t),I{t)) ■±{t) = VR{x{t),I{t)) {-D^u{t,x{t))~^S/R{x{t),I{t)) > 0. 

Combining the above lines we obtain 

■^R{x{t),I{t))i{t) <0, foriG[0,T], 
and hence, thanks to ([SP we deduce 


i(?) > 0, for t G [0,t]. 

Therefore, I{t) is increasing with respect to t and in particular 

i(i) > io) for t £ [0) ‘t]- 

Consequently, from the first line of ([25)1 and ([SP, we obtain that 

R{x{t),Io) > 0, for t G [0,r]. 

It follows that 

x{t) G Do, for i G [0,t]. 

which concludes the proof. □ 

It is then immediate that the initial data given by (|40p verify (14111 . One can verify that the above 
conditions are the only properties that we have used to prove that <I> is well-defined and a contraction. 
Therefore, one can apply again the Banach fixed point Theorem and deduce that there exists a unique 
solution of (l28P - (f2^ for t G [(i — 1)(I, id]. 

4.3 Proof of Lemma 14.21 

(i) We first prove that 


\\V{h) - V(/2)||l-([0,5]xMO < C\\h - l2\\L^im)6. 
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(42) 




Let vi = V{Ii) and V 2 = V{l 2 ), and r = V 2 — vi. From (l36|) we obtain 


dtr = (Vni + VV 2 ) ■ Vr + R{x, I 2 ) — R{x, R), in [0,5] x 
r(0, x) = 0, for all x G 


(43) 


Note that the above equation has a unique classical solution which can be computed by the method 
of characteristics. The characteristics verify 


7 (t) =-Vui(t, 7 ) - Vn2(L7). (44) 

For any (ti,xi) G [0,5] x there exists a unique characteristic curve 7 which verifies 7 (ti) = xi. 
Moreover, this characteristic curve is defined in [0,ti]. 

The local existence and uniqueness of the characteristic curve 7 is derived from the Cauchy-Lipschitz 
Theorem and the fact that Vui is Lipschitz with respect to x, for i = 1,2. The latter property derives 
from the fact that 2L^ < D'^vt < —2Li, for i = 1, 2. 

To prove that the characteristic curve 7 is defined in [0, fi] we must prove that 7 remains bounded in 
this interval. However, this property follows using (I44t) and the fact that, since vi, V 2 ^ B, 

|Vui( 7 )| < Cilyj + 02 , fori = 1 , 2 . 

We are now ready to prove ([38|) . To this end, we multiply (1441) by Vr(f, 7 ) to obtain 

Vr(f, 7 ) - 7 ( 1 ) = -(Vui(f, 7 ) + Vu 2 (f, 7 )) ■ Vr(f, 7 ). 

Combining this with (j43|) we deduce that 

Vr(f, 7 (f)) • 7 (f) + dtr{t, 7 (f)) = R{x, h) - R{x, h). 


We integrate this between 0 and fi to find 





xj2{r)) - R{x,Ii{T))dT. 


It follows, using dSI), that 


and hence ()42D . 


?’(ii,7(ii))l < ^ 2 !^! - 


(ii) Next, we prove that 

II^^IIl°°([0,(5]xR'^) ^ C!\\Ii — I2 \\l°°{[0,S])^- (45) 

We differentiate ([ 43 l) in the direction e^, for f = 1 , • • • , d, to obtain 

dtr(i) = + Vr;2,(i)) • Vr + (Vui + Vr2) • Vr(7 + i?(q(x,/2) - i 2 (i)(x,/i), 

with the notation /(j) = V/ • Cj. We multiply by r(q, sum over i and divide by jVrj to obtain 

n 

SfjVrj < (Vui + VV2) ■ VjVrj + ^ |Vu^(j) + Vr;2,(i)| |r(7| + |Vii’(x,/2) - Vi?(x,/i)|. 

i=l 
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Since D‘^Vj, for j = 1, 2, is bounded, and using ([9]) we deduce that 

dt\Vr\ < (Vui+Vu 2 )-V|Vr| + C|Vr|+i^3|/2-/i|. 

The characteristic curves corresponding to the above equation verify again (Hill . We multiply ()ii|) by 
V|Vr|(t, 7 (t)) to obtain 

V|Vr|(t, 7 ) • 7 (t) = -(Vui(t, 7 ) + Vu 2 (t, 7 )) • V|Vr|(t, 7 ). 

Combining the above equations we obtain 

^|Vr|(t, 7 (t)) < C'|Vr|(t,7(t)) +ii:3||/2 - Ii\\L--i[o,S])- 

It follows that 

|Vr(t,x)| < - l) ^||l 2 -/i||loo([o, 5 ]). 

Hence (j45ll . modifying the constant C if necessary. 

(iii) Finally we prove that 

ll-^^^^lll,°°([0,(5]xR'i) ^ C\\Ii — /2||l°°([0,(5])<5- (46) 

Note that at every point (t, x) G M"*" x we can write mixed derivatives of the form in terms of 
pure derivatives: 

This implies the existence of ^ on the unit sphere of such that 

2 3 

11-^ ^llL°°([0,(5]xR‘i) < 2 

We differentiate (143 jl twice in the direction of ^ and obtain 
dtr^^ = + Vu 2,55) • Vr + 2 (Vui,^ + Vv2,^) ■ Vr^ + (Vui + VU2) • + R^^{x, h) - R^^ix, R). 

Using the above arguments, the fact that D^Vi and D^vi are bounded, ([9]) and (HSl) we deduce that 
^ C\\Ii — -l2||l,oo([0,5])(5 + C'|k^^||L°°([0,5]xK'i) + (^^1 + ^^ 2 ) • V|r^^| + K^\l2 — /l|. 

Next we use the characteristic curves as previously. We multiply ()44|) by V|r^^(t, 7 (t))| and obtain 
V|r^^(t,7(i))l •7(i) = -(Vui(t,7) + Vu2(t,7)) • V|r^^(t,7)]. 

We combine the above equations to obtain 

< C||4l — ^211 LOO ([0,5]) + C'||rgj||Loo([0,5]xRd)- 

We conclude that 

< Ct{\\l2 — III I LOO ([0,5]) + ||r5^||L°°([0,5]xR‘^))- 
Restricting c{T) (and hence 5) if necessary we get 

C6 

lk??lll,°°([0,<5]xR‘*) < “ 4 i||loo([0,5]), 

hence (146)1 . □ 
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5 An example with quadratic uq and R 

Let us study the (instructive) example of a quadratic equation ([T]). In other words we choose 

f uo(x) =-^Aox-x, 

\ R{x,I) =—^Aix ■ X + b ■ X + Iq — I, 

where Aq and Ai are positive definite matrices and b G The Euler-Lagrange equation for the 
dynamic programming principle writes 


f —7 + 2 ^ 17 = 26, 

I - 7 ( 0 ) + 2 ^ 07 ( 0 ) = 0, 7(t) = X. 


(48) 


In order to translate the equivalent system (1281) . we need V-R(/, x) and D‘^u{t,x). The first quantity 
is easily obtained: 

VR{I, x) = —Aix + 6 . 

We then solve (j48p . differentiate with respect to x and denote by T the differential of 7 - the solution 
of (|48P with respect to x, to obtain (after an elementary but tedious computation) 

r(s) = 

B{s)= Id + e-^^^{V^i + 2Ao)-HV^i-2Ao). 


Notice that the result was expected because 7 is linear function of x, and that B{s) is invertible for 
0 < s <t, and the norms of B and B~^ are bounded uniformly in s and t. And so we have 

DMt,x) = -Aor( 0)2 - + Air(s).r(s)^(is := -C{t), (50) 

and the matrix C{t) is bounded away from 0 and + 00 , uniformly with respect to t. This, by the way, 
is not easily seen on the formula (|5Up : the proof of Theorem 11.11 is still what one should use here. The 
equation for x{t) and I{t) is thus straightforward: 

x(t) = C(t)~^{-Aix{t) + 6), 

I{t) = lo - ^Aix{t) ■x{t) + b ■x{t). 

Thus, existence and uniqueness of (x, I) is straightforward. For large t, examination of (l50P and (14^ 
yields 

D‘^u{t,x) ~t^+oo - 

and 

lim x{t) = A7^6, 

t ^+00 ^ 

lim I{t) = Iq H— ARb-b. 

t —^-|-oo 2 

This is consistent with the known behaviour of I{t) and x{t), as well as a closed form of the competition 
increase I{t) — Iq. 



16 





Acknowledgements 


The authors thank G. Barles and B. Perthame for interesting comments on an earlier version of this 
paper, which had a great influence on the present form of this work. They also thank M. Bardi for 
pointing out the possible tractability of the quadratic case. S. Mirrahimi was partially funded by the 
ANR projects KIBORD ANR-13-BS01-0004 and MODEVOL ANR-13-JS01-0009. J.-M. Roquejoffre 
was supported by the European Union’s Seventh Eramework Programme (FP/2007-2013) / ERG 
Grant Agreement n. 321186 - ReaDi - “Reaction-Diffusion Equations, Propagation and Modelling” 
held by Henri Berestycki, as well as the ANR project NONLOCAL ANR-14-CE25-0013. Both authors 
thank the Labex CIMI for a PDE-probability winter quarter in Toulouse, which provided a stimulating 
scientific environment to this project. 

References 

[1] G. Barles, Solutions de viscosite des quations de Hamilton-Jacobi, Mathematiques et Applications, Springer, 
1994. 

[2] G. Barles, L. C. Evans, and P. E. Souganidis, Wavefront propagation for reaction-diffusion systems of PDE, 
Duke Math. J., 61(3) (1990) pp. 835-858. 

[3] G. Barles, S. Mirrahimi, B. Perthame, Concentration in Lotka-Volterra parabolic or integral equations: a 
general convergence result. Methods Appl. Anal., 16(3) (2009) pp.321-340. 

[ 4 ] N. Champagnat, Mathematical study of stochastic models of evolution belonging to the ecological theory of adaptive 
dynamics, Ph.D. Thesis, University of Nanterre (Paris 10), 2004. 

[5] N. Ghampagnat, R. Ferriere, S. Meleard, Unifying evolutionary dynamics: From individual stochastic pro¬ 
cesses to macroscopic models, Theoretical Population Biology, 69(3) (2006) pp. 297-321. 

[6] N. Ghampagnat, R. Ferriere, S. Meleard, Individual-based probabilistic models of adaptive evolution and 
various scaling approximations, 59 (2008) Progress in Probability, Birkhaiiser. 

[7] N. Ghampagnat, P.-E. Jabin, The evolutionary limit for models of populations interacting competitively via 
several resources, Journal of Differential Equations, 261 (2011) pp. 179-195. 

[8] O. Diekmann, Beginner’s guide to adaptive dynamics, Banach Center Publications 63 (2004) pp. 47-86. 

[9] O. Diekmann, P.-E. Jabin, S. Mischler, B. Perthame, The dynamics of adaptation: an illuminating example 
and a Hamilton-Jacobi approach, Th. Pop. Biol., 67(4) (2005) pp. 257-271. 

[10] L. C. Evans and P. E. Souganidis, A PDE approach to geometric optics for certain semilinear parabolic equa¬ 
tions, Indiana Univ. Math. J., 38(1) (1989) pp. 141-172. 

[11] M. Freidlin, Functional integration and partial differential equations, Annals of Mathematics Studies, 109 (1985) 
Princeton University Press, Princeton, NJ. 

[12] M. Freidlin, Limit theorems for large deviations and reaction-diffusion equations, The Annals of Probability, 
13(3) (1985) pp. 639-675. 

[13] Y. Fujita, H. Ishii, P. Loreti, Asymptotic solutions of viscous Hamilton-Jacobi equations with Ornstein- Uh- 
lenbeck operator. Comm. Partial Differential Equations, 31 (2006) pp. 827-848. 

[14] Y. Fujita, H. Ishii, P. Loreti, Asymptotic solutions of Hamilton-Jacobi equations in Euclidean n space, Indiana 
Univ. Math. J., 55 (2006) pp. 1671-1700. 

[15] S. A. H. Geritz, E. Kisdi, G. Meszena, J. A. J. Metz., Evolutionary singular strategies and the adaptive 
growth and branching of the evolutionary tree. Evolutionary Ecology, 12 (1998) pp. 35-57. 

[16] A. Lorz, S. Mirrahimi, B Perthame, Dirac mass dynamics in a multidimensional nonlocal parabolic equation. 
Communications in Partial Differential Equations, 36 (2011) pp.1071-1098. 

[17] S. Mirrahimi, Phenomenes de concentration dans certaines EDPs issues de la biologie, PhD thesis, Universite 
Pierre et Marie Curie, 2011, http://www.math.univ-toulouse.fr/ smirrahi/manuscrit.pdf. 

[18] S. Mirrahimi, B. Perthame, Asymptotic analysis of a selection model with space. To appear in J. Math. Pures 
App.. 


17 



[19] S. Mirrahimi, J.-M. Roquejoffre, Uniqueness in a class of Hamilton-Jaeobi equations with eonstraints, C. R. 
Acad. Sci. Paris, Ser. I (2015), http://dx.doi.Org/10.1016/j.crma.2015.03.005. 

[20] S. Mirrahimi, J.-M. Roquejoffre, Approximation of solutions of selection-mutation models and error estimates, 
in preparation. 

[21] B. Perthame, G. Barles, Dirae eoncentrations in Lotka-Volterra parabolic PDEs, Indiana Univ. Math. J. 57 
(2008) pp. 3275-3301. 

[22] G. Raoul, Etude qualitative et numerique d’equations aux derivees partielles issues des sciences de la nature, PhD 
thesis, ENS Cachan, 2009. 


18 



